Persistence Length Changes DramaticaUy as RNA Folds 
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We determine the persistence length, Ip, for a bacterial group I ribozyme as a function of con- 
centration of monovalent and divalent cations by fitting the distance distribution functions P{r) 
obtained from small angle X-ray scattering intensity data to the asymptotic form of the calculated 
PwLc{r) for a worm-like chain (WLC). The Ip values change dramatically over a narrow range of 
Mg^"*" concentration from ~21 A in the unfolded state (U) to ~10 A in the compact (Ic) and 
native states. Variations in Ip with increasing Na^ concentration are more gradual. In accord with 
the predictions of polyelectrolyte theory we find Ip oc where k is the inverse Debye-screening 
length. 



Elucidating the mechanisms by which RNA 
molecules self-assemble to form three dimensional 
structures is a challenging problem 0,0,13,0]. Because 
the native state (N) cannot form without significantly 
neutralizing the negative charge on 0, |^ the phos- 
phate group, RNA folding is sensitive to the valence, 
size, and shape of the counterions. At low counterion 
concentrations (C) RNA is unfolded (U) in the sense 
that it contains isolated stretches of base-paired stem 
loops that have large dynamical fluctuations. When 
C > Cm, the midpoint of the transition from U to the 
N, RNA becomes compact as a result of formation of 
tertiary contacts. For many RNA molecules, such as 
the Tetrahymena ribozyme and RNase P, folding to the 
native state is preceded by the formation of multiple 
metastable kinetic intermediates (I) 0, S ■ 

The large dynamic conformational fluctuations in 
the U and I states make it difficult to characterize 
their structures. However, small angle scattering ex- 
periments can be used to determine the shape of RNA 
as it folds. The conformation of RNA in the U, N, 
and the I states is characterized by Rg, the radius of 
gyration, and Ip, the persistence length. Small Angle 
X-ray Scattering (SAXS) UlllilllllJ and Small An- 
gle Neutron Scattering (SANS) [l^l experiments have 
been used to obtain Rg as a function of counterions for 
a number of RNA molecules. In contrast, Ip, which is 
a function of C and valence and shape of counterions, 
is more difficult to obtain. 

In this letter, we use SAXS data and theoretical re- 
sults for the worm-like chain (WLC) to obtain Ip for a 
195 nucleotide group I ribozyme from pre-tRNA(Ile) of 
the Azoarcus bacterium as a function of C for monova- 



lent and divalent counterions. The major conclusions 
of the present study are: (i) The experimentally de- 
termined distance distribution functions P{r) can be 
accurately fit using the theoretical results for worm- 
like chains ior r/Rg > 1 where Rg is the radius of 
gyration of RNA. The Ip values, which were calculated 
by fitting P(r) to PwLc{r) for r > Rg, change dra- 
matically from ~ 21 A in the U state to Ip ~ 10 
A in the compact conformation, (ii) The large reduc- 
tion in Ip occurs abruptly over a narrow concentration 
range in Mg^+ whereas the decrease of Ip in Na+ is 
gradual. This result suggests that the compaction of 
RNA resembles a first order transition in the presence 
of multivalent counterions. (iii) For both Na+ and 
Mg^+, the persistence length scales as Ip ~ Ijj where 
Id is the Debye-screening length. From this finding, 
which is in accord with the predictions of polyelec- 
trolyte theory, we find that the intrinsic persistence 
length of RNA is /° ~ 10 A . 

The Azoarcus ribozyine was transcribed in vitro as 
described previously 'l3|- We carried out SAXS mea- 
surements at Argonne National Lab Advanced Photon 
Source (BIOCAT) beamline using 1.05 A X -rays that 
corresponds to ll.SkeV in energy. A sample to detec- 
tor distance of 1.89m allowed us to probe momentum 
transfer (Q) in the range from (-0.007 to -0.266)A"^ 
A quartz capillary flow cell was used to minimize the 
radiation damage due to X-ray exposure of a given 
RNA chain 9]. The measurements at various flow 
rates showed that X-ray radiation damage is negligi- 
ble. Each measurement was averaged from four sepa- 
rate exposures of two seconds each. The SAXS profiles 
were corrected for the background signal which was 
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FIG. 1: (A) Scattering intensity I{Q) as a function of Q for 
195 nucleotide Azoarcus ribozyme at different as a function 
of Mg^"'' concentration in 20mM Tris buffer, pH — 7.5 at 
32°C The values of Mg^"*" concentrations are given in the 
inset. (B) The dependance of on Mg^"*" (squares) and 
Na'^ (triangles) concentration. Solid line is the fit using 
Hill equation in the form A{l-[Mg'^+]'' /{CJ;i, + [Mg'^+]")) + 
Ho, where A,n and yo are adjustable parameters. We find 
n is 3.33 and 1.20 for Mg^+ and Na"*", respectively. 



measured at same buffer concentrations in the absence 
of RNA. Using the background corrected SAXS inten- 
sity as a function of Q, {I{Q)) (Fig. (lA)) the distance 
distribution function, P(r), was calculated by an indi- 
rect inversion algorithm |l4|. The square of the radius 
of gyration is given by i?^ = J r'^P{r)Ar/ J 2P{r)dr. 

Single molecule measurements of RNA subject to 
tension and our analysis of PDB structure have shown 
that the force-extension curves can be fit using a worm- 
like chain (WLC) model Based on these studies 
we assume that RNA is a WLC for which P{r) cannot 
be calculated analytically. However, a simple theoreti- 
cal expression has been derived [3| for the end-to-end 
Re distribution using a mean- field model of WLC. We 
expect that asymptotically (r >> Rg) the behavior of 
P{r) and the distribution of Rg or Re should have 
the same functional form. Thus, for large r we predict 



P(r) should decay as 



PwLc{r) ^ exp 




where x — Ipr/R^. 

At all concentrations of Na+ and Mg^^, Eq. (3) 
fits the data extremely accurately as long as r/Rg > 1 
(Fig. (2)). The excellent fits in Fig. (2) allows us 
to determine Ip as a function of the counterion con- 
centration. When RNA is unfolded at low Na+ or 
Mg2+ concentration, ~ 21 A with i?g ~ 65 A . 
As the concentration of Na"*" increases from about (20 
- 200) mM, Ip gradually decreases. There is a sharp 
decrease in Rg when [Na"*"] ~ 250 mM that is accom- 
panied by a large reduction in the persistence length 
to Ip ~ 10 A. The changes in Ip are even more dra- 
matic in Mg2+ (Fig. (2B)). As Mg^^ increases from 
0.01 mM to 0.26 mM the persistence length changes 
only by about 3 A from Ip ~21 A (0.01 mM) to 
Ip ~ 18.3 A (0.26 mM). In this concentration range 
Rg decreases from 65 A to 60 A. A further increase in 
Mg2+ to 4.26 mM leads to a reduction in Rg to about 
31 A with a dramatic decrease in Ip to about 10 A. The 
near discontinuous change in Rg in Mg^"*" (Fig. (IB)) 
suggests a first-order coil-globule transition in Mg^^. 
While less common in neutral homopolymers, a dis- 
continuous coil-globule transition has been predicted 
to occur in strongly charged polyelectrolytes [T^ . 

To complement the experimental studies we calcu- 
lated the P{r) functions for the native three dimen- 
sional structure of Azoarcus ribozyme using the coor- 
dinates from X-ray crystallography crystal structure 
[l3] (PDB id: 1U6B) and a model based on sequence 
comparison [l^ . The computations were done using 
the coordinates of only the heavy atoms (C, O, P, and 
N) . To compare the results obtained from crystal struc- 
tures and SAXS data, we used only the heavy atom 
coordinates for chain B (excluding nucleotides 1 and 
197) from 1U6B structure to compute P{r). Similarly, 
the exon fragments were excluded from the Westhof 
model. 

The P(r) function from the SAXS data for the 
N state and those obtained using the X-ray struc- 
ture and the Westhof model are in good agreement 
with each other and the SAXS data(Fig. 3A). The 
radii of gyration for the native state calculated using 
= 2^ Ej in - rjf for the X-ray structure 
and the Westhof model are 31.1 A and 30.7 A respec- 
tively. These values agree well with the results from 
the SAXS data {R^ = 30.9 A). The Ip for the native 
state obtained by fitting the crystal structure P{r) to 
Eq. (3) is 11 A, while for the Westhof model we obtain 
Ip ~ 10.8 A. The good agreement between the crystal 




FIG. 2: Distance distribution functions. (A) P{r) func- 
tions at various Mg^^ concentrations at 32° C are obtained 
using Eq. (1) with Q™,i„ = 0.07 A and Q^ax = 0.1 A. The 
solid lines are fits of P{r) to Eq.(3) (B) Same as (A) except 
the counterion is Na^. The concentrations of counterions 
are given in the insets. 



structure (or the Westhof model) and the SAXS mea- 
surement for P(r) and Ip in the Ic state suggests that 
the effects of complexation and interparticle interac- 
tions in the SAXS experiments are neghgible. 

Since RNA appears to be a charged worm-like poly- 
electrolyte, it is of interest to ascertain if the depen- 
dance of Ip on the Debye-screenin g le ngth conforms 
to the theoretical predictions [l^ l2Clj|. The depen- 
dence of Ip on the square of the Debye-screening length 
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S'kIbI where Ib = e /47refcsT is the Bjer- 



rum length and / is the ionic strength) is linear for 
both Na+ and Mg2+ (Fig. (3B)). For both flexi- 
ble [H and stiff polyelectrolytes HI E^l it has 



been shown that Ir, 



JO , jet 

'■p ^ '■p 



where 1° is the intrinsic 



persistence length and the electrostatic contribution is 
Ip^ oc 1/k^. Deviation from the OSF predictions can 
occur for finite-sized flexible polyelectrolytes. How- 
ever, we do not expect such deviations because RNA 
is intrinsically stiff. Surprisingly, over the range of 
Na"*" and Mg^+ concentrations in which the Azoarcus 
ribozyme undergoes the U ^ Iq transition, the exper- 
imental data confirm the predictions of polyelectrolyte 
theory. From the linear fits of Ip to (Fig- (3B)) 
we obtain Z° ~ 10 A which is similar to those found 




FIG. 3: (A) Calculation of P{r) using the coordinates of 
the heavy atoms from the chain B of the crystal structure 
(1U6B), the Westhof model [13], and SAXS data. The 
symbols are given in the inset. The solid lines show the 
fits using PwLcir) in the range 30 A to 70 A . In this 
range the root mean square deviation of the P{r) for the 
1U6B from PwLc{r) is 0.13 A"^ and for the Westhof model 
13] it is 0.10 A~^. The correlation coefficient of the fits 
of P{r) to PwLcir) is 0.98. (B) Dependence of Ip on 1/ 
in Mg^+ (solid circles) and in Na+ (solid squares). 
Lines represent fits to the data. Note that the I/k^ scale 
for Mg^"'' is given on top. 



for single stranded DNA |23ll24j . 

To assess if P{r) for WLC can be used to fit scat- 
tering measurements on other RNA molecules we used 
Eq. (3) and SAXS data for RNase P 9] as a function of 
Mg^+ concentration. Unlike the Azoarcus ribozyme, 
folding of RNase P is best described using three states, 
namely, U, an intermediate I, and the native state, N 
^]. The I state is populated in the Mg^+ range 0.02 < 
Mg2+< 0.2 . From the accurate fit of the SAXS data 
using Eq. (3) for r/Rg > 1, the Ip values are found 
to be 24.5A, 14.lA, and 11. 6A in the U, I, and N 
states respectively. The largest decrease in Ip and the 
associated Rg occurs in the U ^ I transition, which 
is consistent with the notion that the early event in 
RNA collapse is initiated by counterion condensation 
i- 
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FIG. 4: P(r) for RNase P from reference 0. The dis- 
tance distribution function was calculated from I{q) us- 
ing Eq. (1) with a different numerical Fourier Transform 
method|25l| . The U, I and N states are for Mg^^ concen- 
trations 0, 0.1 , and 10 mM respectively. The lines are fits 
using Eq. (3). 

The present work shows that the size and flexibility 
of RNA molecules as a function of counterion concen- 
tration can be obtained using scattering experiments 
and the WLC model. Given that RNA is a highly 
branched and charged polymer, it is surprising that 
the distance distribution functions can be described 
using elasticity-based polymer models for r/Rg > 1. 
Although the structural basis for such a behavior is 
not obvious, the demonstration that single stranded 
DNA 13, double stranded DNA and polypep- 
tide chains [13 also behave like WLC suggests that 
for compatible interactions between biomolecules the 
local flexibility should be similar. 
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